Tight binding parameters for graphene 



Rupali KundrB 

Institute of Physics, Bhubaneswar 751005 India. 

(Dated: July 24, 2009) 

In this article we have reproduced the tight binding tt band dispersion of graphene including 
upto third nearest neighbours and also calculated the partial density of states (due to n band only) 
within the same model. The aim was to find out a set of parameters descending in order as distance 
towards third nearest neighbour increases compared to that of first and second nearest neighbours 
with respect to an atom at the origin. Here we have talked about two such sets of parameters by 
comparing the results with first principle band structure calculation^. 

I. INTRODUCTION 

Graphene is a single sheet of atomic thickness with carbon atoms arranged hexagonally. Though it was an ideal 
two dimensional material of theoretical interest and one of the earliest material on which tight binding band structure 
calculation was done2^, it has triggered recently a lot of interest among people including reinvestigation of many 
earlier results since its experimental discovery in 2004^. Particularly, a large no of people have recalculated the tight 
binding band with nearest neighbour hopping but without overlap integral correctioni^i^&i&ii, some have calculated 
the same by taking into account the overlap integral correction^, out of these only few calculations are there which 
take care of second and third nearest neighbours along with overlap integral corrections^. It is noticed that the first 
nearest neighbour hopping integral (70) lies around 2.5eV-3.0eV when tight binding band is fitted with first principle 
calculation or experimental datai£i£ near the K point of the brillouin zone of graphene but interestingly, when one 
tries to have a good matching of the tight binding band over the whole brillouin zone by including upto third nearest 
neighbour hoppings and overlap integrals, the tight binding parameters arc considered as merely fitting parameters, 
not as physical entities^ i.e, the values of parameters do not decrease consistently as one moves towards second and 
third nearest neighbours. In present work we have fitted our tight binding band with first principle data with an 
objective to get a set of parameters which is free from the above discrepancy and found out a set which gives good 
matching with the first principle data over the whole brillouin zone. Here we have calculated the density of states 
also with second and third nearest neighbours. The article has been arranged in the following manner: section I talks 
about the geometrical structure of graphene; the electronic structure and density of states (DOS) of graphene with 
nearest, next nearest and next to next nearest neighbour couplings arc discussed in section II and finally, section III 
summarizes and concludes over the previous sections. 

II. GEOMETRICAL STRUCTURE OF GRAPHENE 

Since the geometrical structure of a material plays a crucial role in determining the electronic dispersion of the 
material, it is important to look at the details of the structure of graphene before going into the discussion of band 
structure of graphene. The structure of an ideal graphene sheet is a regular hexagonal arrangement of carbon atoms 
in two dimension as shown in fig. 1. It consists of two inequivalent (with respect to orientations of bonds) triangular 
sublattices called A-sublattice and B-sublattice. The unit cell contains one A and one B type of carbon atoms con- 
tributed by respective sublattices. Each carbon atom has three nearest neighbours coming from the other sublattice, 
six next nearest neighbours from the same sublattice and three next to next nearest neighbours from the other sublat- 
tice. ao (1.42A) is the nearest neighbour lattice distance. In the figure, a\ and 0,2 are the unit vectors with magnitude 
a = V3ao i.e. 2.46A. With respect to Aa atom the coordinates of the first neighbours (Bu,i — 1,2,3); second 
neighbours (A-M,i — 1, ...6) and third neighbours (i?3i,i = 1,2,3) are (a/V3, 0) , (-a/2\/3, -a/2) , (-a/2\/3, a/2) ; 
(0, a) , (V3a/2, a/2) , (VZa/2, -a/2) , (0, -a) , (-\/3a/2, -a/2) , (- V3a/2, a/2) and (a/VS, a) , (a/V3, -a) , 
(— 2a/-\/3, 0) respectively. 



III. ELECTRONIC STRUCTURE OF GRAPHENE 



Carbon atom has six electrons with the electronic configuration ls 2 2s 2 2p 2 . 2s and 2p levels of carbon atoms can mix 
up with each other and give rise to various hybridized orbitals depending on the proportionality of s and p orbitals. 



FIG. 1: Structure of graphene. The structure with black circles forms the A-sublattice and that with red circles gives B- 
sublattice. a± and S2 are the unit vectors. 



Graphene has sp 2 hybridization: 2s orbital overlaps with 2p x and 2p y orbitals and generates three new inplane sp 2 
orbitals each having one electron. The 2p z orbital remains unaltered and becomes singly occupied. Due to overlap 
of sp 2 orbitals of adjacent carbon atoms strong a (bonding) and a* (antibonding) bonds are formed. The bonding 
cr bonds, lying in a plane, make an angle of 120° with each other and is at the root of hexagonal planar structure 
of graphene. p z orbitals being perpendicular to the plane overlap in a sidcwisc fashion and give ir (bonding) and tt* 
(antibonding) bonds, sp 2 orbitals with a lower binding energy compared to Is (core level) are designated as semi 
core levels and p z orbitals having lowest binding energy are the valence levels. Overlapping of p z energy levels gives 
the valence band (bonding 7r band) and conduction band (antibonding tt* band) in graphene. Thus, we see that 
while the structure of graphene owes to <j bonds, 7r band is responsible for the electronic properties of graphene and 
hence as far as electronic properties of graphene are concerned, concentration is given only on tt bands. Since the p z 
orbitals overlap in a sidewise manner, the corresponding coupling is weaker compared to that of a bonds (where sp 2 
orbitals overlap face to face). So the p z orbitals almost retain their atomic character. Hence, to describe the electronic 
structure of graphene, tight binding model could be a good choice. 



A. Tight binding band 

In this subsection tight binding bands of graphene have been reproduced including upto third nearest neighbor 
hopping of electrons and overlap integral corrections with focus on the point to find out tight binding parameters 
which will not look like just fitting parameters rather look somewhat like physical parameters. The results have been 
compared with existing literature^. The key equations are shown here for the different cases of nearest neighbours, 
next nearest neighbours and next to next nearest neighbours and the details of the calculations are given in the 
appendix. Since there are two atoms per unit cell coming from two sublattices, the total wave function can be written 
as 

f fc (r) = C A *\ + C B *%, (1) 

where ^ A (r) = l/\/N'E A e ik - r ~ A <& A (r - r A ) and &%(r) = l/VNE B e ikr ~ B § B (r - Tb) arc the tight binding Bloch wave 
functions from A and B sublattices. Here, N is the number of unit cells in the crystal, C A and C B are contributions 
coming from A and B sublattices respectively, <E> s are 2p z atomic orbitals, k is crystal momentum and r A and r B are 
the positions of A and B atoms respectively with respect to a chosen origin. If H is the Hamiltonian and E(k) the 
eigenvalue then 



(2) 
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which leads to the secular equation 



H AA - E{k)S AA H AB ~ E{k)S AB 
H* AB - E(k)S* AB H BB - E(k)S AA 







Since two sublattices are equivalent, H AA = Hbb and S AA = Sbb and the general dispersion relation follows as 
E ± (k) = \(2E - Eh.) ± V(Ei -2E ) 2 - ±E Y E 2 



(3) 



(4) 



>AA 



SabS 



AB 



E 3 . The explicit 



where S AA H AA = E , H AB S* AB + H* AB S AB = E u H 2 AA - H AB H* AB = E 2 and S'j 
forms of the dispersions under various situations are discussed below. 

With Nearest Neighbour Approximation: In case of nearest neighbour approximation contribution comes from 
nearest atoms of the other sublattice. Detailed calculation of all the matrix elements is shown in appendix. From 
equation (4) the dispersion relation for this case becomes 



E ± (k) = (E 2p - s 7o<?(fc)) ± (7o - s E 2p ) ^/gjk) / [l - s 2 g(k)] , 



(5) 



where E 2 p, 70 and sq are site energy, nearest neighbour hopping and overlap integrals respectively. 
With Second Nearest Neighbour Approximation: In this part the modification in the energy spectrum due to 
the presence of second nearest neighbor atoms will be discussed. The matrix elements H AA and S AA get modified 
but H A b and Sab remain as they were. The corresponding energy momentum relation is 



E ± (k) 



j 2p 



7iit(fc) T 7o\/sW / 1 + siu(k) =p s oy / g(k) 



(6) 



where 71 and s± are next nearest neighbour hopping and overlap integrals respectively. 

With Third Nearest Neighbor Approximation: In this case the matrix elements H A b and Sab get changed but 
Haa and S AA remain unaltered. Calculation of all the matrix elements gives expressions for E\, E 2 , E3 as follows: 

Ei = 2s 7off(fc) + (s 72 + 70S2) t(k) + 2s 2 j 2 g(2k) 
E 2 = [E 2p + 7 iw(fc)] 2 - [i 2 g{k) + 7 o72t(fc) + i 2 2 g{2k)} 
E 3 = [1 + Sl u(k)} 2 - [s 2 Q g(k) + s s 2 t(k) + s 2 g{2k)} , 

which when put in equation (4) gives the energy momentum relation in this case. 72 and s 2 being next to next 
nearest neighbour hopping and overlap integrals respectively. The bands of graphene for above three cases are 
plotted along the high symmetry directions of its hexagonal brillouin zone sketched in figure (2). As shown in 




FIG. 2: Brillouin zone of graphene. (a) Real space unit vectors, (b) k-space unit vectors and high symmetry directions. 



the figure, the symmetry points are r(0, 0), M(2vr/\/3a, 0) and K(2n/y/3a,2Tv/3a). There are six corner points, out 
of which three are independent since the nearby corner points(-ftT and K') are inequivalcnt. In fig. (3) we have 
compared the bands with nearest, next nearest and next to next nearest neighbour hopping and overlaps with a first 
principle calculation [produced from this model with the parameters E 2 p = — 0.36cV, 70 = — 2.78cV, 71 = — 0.12cV, 
72 = — 0.068eV, so = 0.106, si = 0.001, s 2 — 0.003 because within the chosen energy scale it does not show any 
energy difference with first principle datai] for the set of parameters shown in table I where we have first determined 
values of nearest neighbour parameters (70 and sq) which best reproduces the first principle result, then we have 



r km r 
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FIG. 3: Electronic structure of graphene with first principle result (black curve), first nearest neighbour interactions (green 
curve), second nearest neighbour interactions (blue curve) when first nearest neighbour parameters are fixed and third nearest 
neighbour interactions (red curve) when first and second nearest neighbour parameters are fixed. The parameters are listed in 
Table I. 



TABLE I: Tight binding parameters 



Neighbours 


E 2P (eV) 


70 (eV) 


7i (eV) 


72 (eV) 


so 


si 


S2 


fst 


0.0 


-2.74 






0.065 






2nd 


-0.21 


-2.74 


-0.07 




0.065 


0.002 




3rd 


-0.21 


-2.74 


-0.07 


-0.015 


0.065 


0.002 


0.001 



taken care of second nearest neighbours (71 and s\) keeping first nearest neighbour parameters fixed and lastly 
considered the third nearest neighbour parameters (72 and S2) with fixed (70, so) and (71, si) with the expectation 
to have better matching over the whole brillouin zone of graphene. We observe that nearest neighbour coupling 
gives good overall matching for both valence and conduction bands but with the inclusion of second and third 
nearest neighbours under the above restrictions, the overall matching of the valence band with the first principle 
band improves whereas for the conduction band it is good in the optical range. When only nearest neighbour 
hopping is considered, the total band width (difference between valence band and conduction band energies) at T 
point is 6|7o|, that at M point is 2|7o| but when nearest neighbour overlap integral is included the valence and 
conduction band energies at T point occur at 370/(1 + 3so) and — 370/(1 — 3so) respectively and those at M point 
appear at 7o/(l + so) and —70/(1 — so) respectively. The K point energy is zero for both the cases. When next 
nearest neighbours are included the energy at K point is (E 2p — 37i)/(l — 3si). So the values of E 2p and 71 
are properly chosen to have K point energy at zero. In this case the valence and conduction band energies at T 
point are at {E 2p + 6si + 37o)/(l + 6si + 3so) and (E 2p + 6si — 37o)/(l + 6si — 3so) and those at M point are at 
(E 2p + 2si + 7o)/(l + 2si + So) and (E 2p + 2si — 7o)/(l + 2si — s ) respectively. In fig. (4) we have plotted the 
above set of curves with a different set of parameters such that both nearest neighbour and next nearest neighbour 
parameters are free when the effect of second nearest neighbour is looked for and all three, nearest, next nearest and 
next to next nearest neighbour parameters are free when the third nearest neighbours are included. The parameters 
are shown in table II. This figure shows that inclusion of second nearest neighbour gives better result over first nearest 
neighbour interaction compared to the plots in fig. 3. Consideration of third nearest neighbour couplings gives very 
good matching for both valence and conduction bands along all the high symmetry directions of the brillouin zone. 



TABLE II: Tight binding parameters 



Neighbours 


£ 2 p(eV) 


70 (eV) 


71 (eV) 


72 (eV) 


so 


Sl 


S2 


1st 


0.0 


-2.74 






0.065 






2nd 


-0.30 


-2.77 


-0.10 




0.095 


0.003 




3rd 


-0.45 


-2.78 


-0.15 


-0.095 


0.117 


0.004 


0.002 
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FIG. 4: Electronic structure of graphene with first principle result (black curve), first nearest neighbour interactions (green 
curve), second nearest neighbour interactions (blue curve) and third nearest neighbour interactions (red curve). The parameters 
are listed in Table II. Here for each curve the parameters are chosen freely. 



B. Density of states of graphene 

The density of states of graphene for the above two cases are shown in figures (6) and (7). In fig. (5) we have 



< 




Energy in eV 



FIG. 5: Density of states of graphene with nearest neighbour hopping but zero overlap integral (black curve) and with nearest 
neighbour hopping and overlap integral (red curve). 

plotted the density of states of graphene (only for tt band) for nearest neighbour interaction only without (black curve) 
and with (red curve) overlap integral correction. It shows that without overlap valence and conduction bands are 
symmetric but in presence of overlap width of valence band decreases and that of conduction band increases. Since 
there are two atoms in each unit cell and each carbon atom has one electron in p z state, the valence band is completely 
filled and hence, the fermi level lies at the top of the valence band at zero energy which appears at K and K' points 
of the brillouin zone in energy momentum space. Graphene has zero density of states at fermi energy and over a very 
small energy range around zero (i.e. K point energy around which energy dispersion is also linear in momentum) the 
density of states are varying linearly with energy. Due to the flat part of the band near M point of brillouin zone 
van-Hove singularities are arising in density of states. The positions of the singular points are symmetric in absence of 
overlap term while the singularity moves slightly towards fermi energy for valence band and goes slightly away from 
fermi energy for conduction bands with overlap (sq). In figure (6) and (7) the density of states have been plotted 
as per the bands in fig. (3) and (4). Here we see that though the density of states is linear near fermi energy, it 
becomes asymmetric due to the presence of second and third nearest neighbour interactions. As expected the density 
of states in presence of third nearest neighbours is matching well with the density of states due first principle bands 
when the parameters are chosen freely. In all other cases the band width is changing slightly and positions of van-Hov 
singularities are changing over a narrow energy range. When compared with rcf . (9) , the density of states in presence 
of second nearest neighbours shows a noticable difference. It is observed that when the second nearest neighbours are 
chosen for good fitting of low energy part of the spectrum, it gives rise to another van-Hov singularity at conduction 
band edge but that is removed if the parameters are such that it gives good fitting over the whole energy range. Also, 
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FIG. 6: Density of states of graphene derived from the bands in fig. 3 (parameters are in Table I). First principle result (black 
curve), first nearest neighbour interactions (green curve), second nearest neighbour interactions (blue curve) and third nearest 
neighbour interactions (red curve). 




FIG. 7: Density of states of graphene calculated from the bands plotted in fig. 4 (parameters are in Table II). First principle 
result (black curve), first nearest neighbour interactions (green curve), second nearest neighbour interactions (blue curve) and 
third nearest neighbour interactions (red curve). 



IV. SUMMARY AND CONCONCLUSION 



To summarise, we have calculated the electronic spectrum and density of states of graphene including upto third 
nearest neighbour interactions and got a set of tight binding parameters on the physical ground that the absolute 
values of the parameters should decrease as one moves from first nearest neighbour towards higher distance. This set 
of parameters can be used to see the effect of intraplane first nearest neighbour overlap integral (so) and second (71, 
si) and third (72, S2) nearest neighbour interactions on the band structure of bilayer graphene and graphite. 
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APPENDIX A: DERIVATION OF THE ELECTRONIC DISPERSIONS IN GENERAL AND IN 

PRESENCE OF DIFFERENT NEIGHBOURS 



In this appendix we present slightly more details of the derivations of dispersion relations in presence of various 
nearest neighbours. We have the wave function 



ttfc(r) = C A * k A + C B ^ k B 
If H is the Hamiltonian and E(k) the eigenvalue then 

HV k (r) = E(k)V k (r) 
which gives the secular equation 



H AA - E(k)S AA B AB - E(k)S AB 
H* AB - E(k)S* AB H BB - E(k)S AA 

where H xy = J ^f{r)H^ h y {r)dr 







(Al) 
(A2) 

(A3) 



H 



and Sxy — 



S u 



where x, y represent A and B both. Since two sublattices in graphene are equivalent El a A = Hbb and Saa = Sbb- 
.'. the general dispersion relation is 



E ± (k) = (2E a - E{) ± ^{E 1 - 2E ) 2 - ±E X E 2 /2E 3 , 



(A4) 



where we define SaaHaa = E , HabS ab + H ab S A b = E x , H AA - H AB H AB = E 2 and S\ A - S AB S* AB = E 3 . 
With Nearest Neighbour Approximation: Here we calculate the elements of the determinant in equation (A3) 
for the nearest neighbour interactions of the atoms. 



Haa 
Saa 

Hab 

where 70 
/(*) 



l/NSAe^-^ ($ A (r - r A )\H\$ A (r ~ r A )) « E 2p . 

(^a(t — r A )\$> A (r — r A )) = 1 (since the wave functions are normalized). 

7o/(fc), 

1/NTiA (^a(t — ta)\H |$s(r — rg)) is the hopping integral (-ve) and 

e ik„a/V5 + 2e -ik*a/2V5 CQS kyCL /2, 



Similarly, Sab 
where sq 
:.E 
Ei 
E 2 
and E3 
Also, |/(fc)| 2 



sof(k), 

\/NYja (3>a(»" — ta)^b(i" — i'b)) is the overlap integral (+ve). 
E 2p 

27os |/(fc)| 2 
E 2 2p -J 2 o\f(k)\ 2 
l-4\f(k)\ 2 
9(k) 

1 + 4 cos 2 (kyd/2) + 4 cos (V3k x a/2) cos (k y a/2) 
(E 2p - so-fog(k)) ± (70 - s E 2p ) y/g(k) / [l - slg(k)] . 



An Alternative Expression: From the secular determinant the dispersion can alternatively written as£ 



E ± (k) = (E 2p ± l0 ^gjk)) I (l ± soVsM) 



(A5) 



(A6) 
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With Second Nearest Neighbour Approximation In this part the matrix elements Haa and Saa are calculated 
because only these two elements are affected in presence of next nearest neighbours. 

H AA = J y*A k (r)HV A {r)dr 

= E 2p + Jiu(k), 

where u(k) = 2 cos (k y a) + 4 cos (k x aV3) cos (k y a/2) and 

71 = 1/NH A > {& A > (r — r A ')\H\&A{r — ta)) is the second neighbour hopping integral. 

Similarly, Saa = 1 + s\u(k) 

si = l/NT, A ' (& A >(r — r A i)\<^ A {r — ta)) is the second neighbour overlap integral. 

Putting all the above elements in the secular determinant we get 



E ± (k) = \E 2p + 7iu(fc) T 7o\/flWj / [l + stu(k) + s ^/gjk)\ . (A7) 

With third nearest neighbor approximation: In presence of third nearest neighbours only the matrix elements 
Hab and Sab get changed but Haa and Saa remain unchanged. Here 

H AB = J K B k (r)HV A (r)dr 

= 7o/0) +l2v(k), 

where72 = 1/NY, A ~ Ta)\H\& b {t — t b)) is the third nearest neighbour hopping integral 

andi»(fc) = e^ a/V ^2cosk y a + e~ 2ik " a ^. 

Similarly, Sab = s f(k) + s 2 v(k), 

where s 2 = l/NT, A {$ A (r - r A )\$B(r - r B )) ■ 

.'. expressions for E±, E 2 , E$ are 

Ei = 2so7o.9(fc) + («o72 + 70S2) t{k) + 2s 2 j 2 g(2k) 
E 2 = [E 2p + Jiu(k)] 2 - [^g(k) + 7 o72i(fc) + llg(2k)} 
E 3 = [1 + s lU (k)} 2 - [s 2 g(k) + s s 2 t(k) + s 2 2 g(2k)] , 
where g(2k) = 1 + 4 cos 2 ik y a) + cos (V3k x a) cos (k y a) 

t(k) — 2 cos (k x aV3) + 4 cos (k y a) + 4 cos (k y a/2) cos (k x aV3/2) 
+8 cos (k y a) cos (k y a/2) cos (k x aV3/2). 
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